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$8 要 : 大 气 水 汽 氨 氧 稳定 同位 素 直接 反映 了 水 在 大 气 中 输送 混合 和 相 变 等 过 程 的 关键 信息 ,通过 蒸发 四 实验 
可 以 模拟 大 气 水 汽 氨 氧 稳定 同位 素 组 成 ,但 模拟 方法 的 可 靠 性 仍 需 要 实测 资料 的 验证 。 于 2019 年 9 一 11 月 在 甘肃 


兰州 进行 了 蒸发 严实 验 , 对 逐日 大 气 水 汽 氧 氧 稳定 同位 素 进行 模拟 ,同期 利用 在 线 大 气 水 汽 氧 氧 稳定 同位 素 分 析 
仪 进行 实时 监测 从 而 对 模拟 结果 的 可 靠 性 进行 验证 。 研 究 发 现 :(1) 利用 自然 蒸发 状态 下 的 Craig-Gordon 线性 阻力 
模型 和 水 体 同 位 素质 量 平衡 方法 模拟 得 到 的 大 气 水 汽 中 3"0 与 直接 测量 得 到 的 大 气 水 汽 中 3*0 之 间 表 现 出 较 好 的 
一 致 性 , 均 方 根 误差 为 4.5%o ,平均 绝对 误差 为 3.2%o ,平均 偏差 为 0.03%o。(2) 大 气 水 汽 中 50 的 模拟 残 差 与 相对 湿 
度 (R=0.43) .蒸发 剩余 比 ( 尼 =0.39) 和 剩余 水 体 同 位 素 值 ( 尼 =0.39) 都 具有 一 定 相关 性 。(3) 模型 对 相对 湿度 的 敏感 
性 高 ,在 相对 湿度 较 低 时 模型 误差 会 增 大 。 在 实验 末期 薰 发 剩余 比 偏 低 和 剩余 水 体 同 位 素 富 集 的 情况 下 ,水 汽 同 


位 素 模拟 残 差 相 较 于 实验 前 期 而 言 更 大 。 
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水 体 中 的 氧 氧 稳定 同位 素 (H 和 ?0) 是 研究 水 
文 循环 和 十 气候 信息 重建 的 天 然 示 踪 剂 。 大 气 
水 汽 是 形成 降水 的 必要 物质 条 件 ,水汽 中 氧 氧 稳定 
同位 素 也 相应 地 列 含 着 蒸发 .输送 .凝结 等 环境 信 
息 “。 相 比 于 降水 ,大 气 水 汽 同 位 素 能 为 水 汽 过 程 
诊断 提供 直接 的 证 据 , 其 优势 表现 在 时 间 空间 等 
多 方面 5。 水 汽 同 位 素 可 以 实现 全 天 候 的 连续 监 
W= ,不 像 降水 同位 素 那样 仅 能 够 代表 降水 时 段 
的 特征 ,从 而 得 到 时 间 序 列 连续 的 大 气 水 分 信息 ， 
这 在 降水 匮乏 的 干旱 半 干 旱 区 显得 十 分 重要 。 结 
合 浮 空 艇 .无 人 机 和 卫星 遥感 等 手段 ,水 汽 同位 素 
还 可 以 实现 三 维 空间 覆盖 , 相 比 之 下 地 面 降 水 同位 
素 对 垂直 方向 上 的 水 汽 过 程 刻 画 则 具有 其 局 限 
性 "中 。 此 外 ,雨滴 在 不 饱和 空气 中 的 下 落 过 程 可 
能 会 发 生 云 下 二 次 蒸发 ,影响 降水 同位 素 的 时 空 特 
PEC ,而 对 水 汽 同 位 素 的 监测 能 够 在 一 定 程 度 上 
弱化 云 下 莹 发 导致 的 降水 同位 素 信 息 失 真 问题 。 

早期 大 气 水 汽 氧 氧 稳定 同位 素 研 究 主要 以 低 
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温 冷凝 空气 获得 液态 水 样 进行 分 析 , 即 低温 冷 阱 
法 ,该 方法 对 实验 操作 具有 较 高 的 要 求 ,往往 难以 
精细 表征 大 气 水 汽 同 位 素 在 时 间 和 空间 上 的 变化 ， 
特别 是 在 水 汽 含量 低 的 情况 下 水 汽 液化 比较 耗 时 
费力 ”"。 随 着 激光 光谱 技术 的 发 展 ,实时 在 线 监 
测 大 气 水 汽 的 氧 氧 稳定 同位 素数 据 成 为 可 能 ”。 
激光 光谱 技术 克服 了 低温 冷 阱 法 在 连续 测量 上 的 
局 限 性 ,但 由 于 成 本 相对 较 高 ,无 法 广泛 地 开展 监 
测 工作 。 当 不 能 直接 监测 大 气 水 汽 同 位 素 时 ,可 利 
用 Craig-Gordon 线 性 阻力 模型 ”和 水 体 同位 素质 量 
平衡 方法 来 模拟 ”1 ,并 结合 实测 水 汽 同位 素数 据 
加 以 验证 。 

为 了 人 研究 基于 Craig-Gordon 线 性 阻力 模型 和 同 
位 素质 量 平衡 的 方法 在 计算 干旱 半 干 旱 区 大 气 水 
汽 毛 氧 稳定 同位 素 时 的 可 靠 性 ,本 文 基于 燕 发 甩 实 
验 模拟 了 甘肃 兰州 大 气 水 汽 氧 氧 稳定 同位 素 ,并 结 
合 在 线 水 汽 同 位 素 分 析 仪 的 同步 监测 记录 ,评估 了 
燕 发 邮 实 验 模拟 结果 的 可 靠 性 , 旨 在 为 获取 连续 大 
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气 水 汽 氢 氧 稳定 同位 素数 据 提供 了 一 种 相对 简单 
的 方法 。 


1 材料 与 方法 


1.1 AAO 
试验 地 位 于 甘肃 省 兰州 市 安宁 区 西北 师范 大 
FIREX (103?44' E, 3626' N) ,地 处 季风 边缘 区 , 属 
于 半 干 早 气 修 , 气 温 年 较 差 大 ,降水 多 集中 在 夏 
Æ, EAMUS F 201949 H25 H—11 H 1 H3 
行 , 实 验 期 间 的 平均 气温 为 14.4 C ,平均 相对 湿度 
为 46%。 选 用 长 420 mm, $ë 220 mm 高 200 mm 的 
长 方 体 透 明 玻 璃 水 箱 作 为 蒸发 四 ,放置 于 开阔 388 
风 良 好 的 场地 。 当 有 降 十 发生 时 ,在 蒸发 阵 上 方 
0.5 mm 高 度 处 用 透 光 塑料 薄膜 遮盖 ,防止 降雨 进入 蒸 
A MIL ,降雨 结束 后 及 时 移 除 。 实 验 初始 将 当地 自 来 
水 注入 蒸发 阵 中 ,每 日 同一 时 间 用 钢 尺 测量 水 位 并 
记录 ,然后 充分 搅拌 水 体 后 用 注射 器 抽取 1 mL 的 水 
体 样 本 储存 在 玻璃 瓶 中 密封 冷藏 保存 。 为 了 减少 
取样 造成 水 量 损失 对 结果 的 影响 ,在 上 述 蒸发 血 旁 
0.5 m 处 设置 一 个 同样 大 小 和 材质 的 燕 发 呈 作 为 补 
给 箱 ,在 初始 条 件 下 装 入 同 量 的 自来水 同时 进行 蒸 
发 ,每 日 抽取 1 mL 的 水 体 补充 到 用 于 收集 剩余 水 体 
的 蒸发 生 中 。 实 验 过 程 中 实时 监测 大 气温 度 和 相 
对 湿度 并 记录 。 抽 取 的 逐日 蒸发 剩余 水 体 采 用 加 
拿 大 ABB-Los Gatos Research 公司 人 研发 的 T-LWIA- 
A5—EP 型 液态 水 同位 素 分 析 仪 进行 氧 氧 稳定 同位 素 
测定 ,5*0 的 测量 误差 小 于 +0.3%o , SH 的 测量 误差 
小 于 +1.0%o。 
氧 氧 稳定 同位 素 比 率 用 相对 于 维也纳 标准 平 
均 海洋 水 (V-SMOW ) 的 千 分 差 表示 ,公式 如 下 : 
D NES Re 
= | sample V-SMOW 


R V-SMOW 


» 1000966 (1) 


TUF : Roane RIR TRACK SEAS HY FE R ELK ; Rv swow 
表示 维也纳 标准 平均 海洋 水 的 同位 素 比 率 。 
12 大 气 水 汽 同 位 素 测量 

大 气 水 汽 氨 氧 稳定 同位 素 的 监测 采用 加 拿 大 
ABB-Los Gatos Research 公司 研发 的 超 便携 式 水 汽 
同位 素 分 析 仪 (Ultraportable Water Vapor Isotope An- 
alyzer, 型 号 UWVIA-915) , 55 Z& Az IML Sz 3 [vi] 80] ff 
行 。 该 分 析 仪 基于 离 轴 集 成 腔 输出 光谱 技术 ,可 以 
通过 测定 不 同 氨 氧 稳定 同位 素 对 于 特定 波长 红外 


波 的 吸收 强度 来 计算 其 浓度 ,采样 频率 可 达 秒 级 。 
本 研究 中 8 "0 的 测量 精度 为 40.05%o, SH 的 测量 精 
度 为 +0.2%o。 仪 器 通过 水 汽 校准 源 系统 (Water Va- 
por Isotope Standard Source, 简称 WVISS) 产 生 标 样 
水 汽 , 测 量 过 程 中 将 标 样 水 汽 每 浓度 梯度 开始 时 和 
结束 时 的 测量 值 删除 , 取 每 梯度 中 间 较 为 稳定 的 同 
位 素 值 作为 该 梯度 的 标 样 水 汽 同 位 素 值 以 减轻 记 
忆 效 应 对 实验 结果 的 影响 。 在 测量 大 气 水 汽 同 位 
素 时 ,将 前 10 min 的 测量 值 删 除 ,至 大 气 水 汽 浓 度 
稳定 以 后 再 进行 测量 »”。 浓 度 效 应 的 校准 参考 了 
Rambo 等 ”的 方法 , 标 样 水 中 同位 素 和 气 化 后 的 标 
样 水 汽 中 的 同位 素 一 致 ™ ,可 通过 T-LWIA-45-EP 
型 液态 水 同位 素 分 析 仪 测定 后 作为 标 样 水 汽 的 真 
值 ,将 该 测量 真 值 与 水 汽 同位 素 分 析 仪 得 到 的 标 样 
水 汽 同位 素 求 差 , 用 水 汽 浓度 和 差 值 进行 二 项 式 曲 
线 拟 合 ,获得 对 应 水 汽 浓度 下 的 测量 误差 值 ,从 而 
对 水 汽 同 位 素 进行 浓度 效应 校准 。 
13 基于 蒸发 血 实 验 的 模拟 计算 

Craig-Gordon 线性 阻力 模型 被 广泛 应 用 在 薰 
发 水 汽 同位 素 的 估算 中 ,方法 如 下 : 
3,-e)/a - h8, 
: 1 zt Ex (2) 
SUF ioi ZR OKITA HY IE LR LOR 5 8, REN ZR HE 
余 水 体 的 同位 素 比 率 ;5, 表 示 大 气 水 汽 同位 素 比 率 ; 
hh 是 空气 的 相对 湿度 ;a 是 平衡 分 馏 系 数 ,可 以 通过 
Horita 等 ? 的 公式 得 到 ;es 是 平衡 分 馏 因 子 , 可 由 s = 
a -1 得 到 ;sx 是 动力 分 馅 因子 , 氧 和 氢 的 同位 素 动力 
4 RES 49] ex=14.2(1-h) Fllex=12.5(1-h) ; & 3b M 
4T BAT esee, Gibson E WHA BIKA 
同位 素 取决 于 大 气 水 汽 的 同位 素 ,对 于 一 个 混合 
好 的 水 体 在 自然 燕 发 条 件 下 而 言 , 在 di 的 时 间 间 隔 
下 ,水 量 的 变化 为 dV, 则 质量 平衡 方程 为 : 


8,- 


X =1-Q-E (3) 
同位 素质 量 平衡 方程 为 : 
309) — p Q5, - pi, (4) 


式 中 :V 是 水 体 的 水 量 ;7 是 水 的 流入 量 ;5 是 流入 水 
体 的 同位 素 组 成 ;0 是 流出 量 ;5 是 流出 水 体 的 同位 
素 比 率 ;E 是 蒸发 量 ;6: 是 蒸发 水 体 的 同位 素 比 率 。 
将 式 (2) 中 8: 代 入 式 (4) 可 得 : 
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(8, - e)/o -hð, -£x 
(l-h+e,) 


yo. " s, 4K -[8,- Qò,- E 


(5) 

此 式 为 同位 素质 量 平衡 方程 的 通用 公式 。 在 

本 实验 中 , 蔡 发 亚 作 为 蒸发 水 体 的 容器 ,体积 固定 ， 

不 存在 流入 和 流出 的 情况 , 落 发 呈 中 水量 只 受到 自 

然 条 件 下 的 藻 发 影响 而 减 小 。 将 WW 当 作 蒸发 初始 

时 的 水 量 ,定义 庆 凡 有 为 蒸发 程度 的 指标 , 式 (4) 可 
改写 为 : 


dò, _ | 
dünf) 9, (6) 
将 式 (2) 和 式 (6) 合 并 整理 可 得 : 
6,-8 — i 
Z (7) 
y Rte (8) 
h-e 


BERERE 9 ET KP KO HOT 1 
个 固定 值 , 称 之 为 同位 素 的 “ 稳 态 ”, 即 式 (7) 中 8 ,其 
中 m= 于 2 一 是 一 个 几乎 完全 取决 于 相对 湿度 的 


将 式 (7) 和 式 (8) 整 理 后 可 得 : 
Bs je -8,f")(h - e) -e (9) 
NER 

TAX BI Jy 3 ER e ILS TT OK TM 
素 的 计算 式 。 

对 于 大 气 水 汽 同位 素 的 模拟 精度 采用 均 方 根 
误差 (Root Mean Square Error, RMSE) .平均 绝对 误 
差 (Mean Absolute Error, MAE ) 和 平均 偏差 (Mean Bi- 
as Error, MBE ) 进 行 评估 。 


(10) 
MAE = Tòni- b (11) 
MBE- 15 61i- Srm.) (12) 


E 中 让 表示 实测 的 大 气 水 汽 同位 素 比率 ; Sm 表示 
模拟 的 大 气 水 汽 同位 素 比 率 。 


2 结果 与 分 析 


2.1 和 氢 氧 稳定 同位 素 的 时 间 变 化 
蒸发 水 体 的 稳定 同位 素 初始 值 分 别 为 8*0= 
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-9.50%o 和 82H=-65.60%o , 随 着 蒸发 阵 中 水 体 的 蒸发 
损失 ,剩余 水 体 中 稳定 同位 素 逐 渐 富 集 ,350 从 
-9.50%o 增 大 至 21.08%o( 图 la) ,5H 从 -65.60%o 增 大 
至 50.25%o( 图 1b)。 这 符合 蒸发 剩余 水 体 中 同位 素 
的 基本 变化 规律 只 35。Devi 等 避 在 印度 的 蒸发 自 实 
验 中 ,剩余 水 体 中 80 从 -6.30%o 增 大 至 13.63%o,5°H 
从 -43.48%o 增 大 至 46.50%o ,与 本 研究 中 水 体 同位 素 
的 变化 范围 比较 接近 ,但 实验 期 间 平 均 温 度 
(26.5 % ) 和 相对 湿度 (49% ) 比 本 研究 更 高 。 

在 实验 前 期 ,蒸发 剩余 水 体 同 位 素 变化 相对 较 


0.8 


EEETINO 
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图 1 20194 9—11 HRE MKS SCRI ZR ACHR AK PI 
(SL) ,模拟 蒸发 水 汽 (6:) .模拟 大 气 水 汽 (546w) 和 实测 大 气 水 
汽 (56,) 中 5"*O(a) 和 8*H(b) 的 日 变化 
Fig. 1 Daily variations of O (a) and 6°H (b) in measured 


remaining pan water (6,), simulated evaporation water vapor 
(05), simulated atmospheric water vapor (04.55), measured 
atmospheric water vapor (õa) during evaporation pan 


experiment from September to November 2019 
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为 平缓 ,而 在 实验 末期 蒸发 血 中 水 量 减 少 ,水 体 更 
易 受 到 气象 条 件 变 化 的 影响 ,同位 素 的 芋 发 富 集 速 
率 略 微 增加 。 伴 随 着 蒸发 过 程 ,模拟 的 蒸发 水 汽 同 
位 素 也 逐渐 富 集 ,表现 出 和 蒸发 剩余 水 体 类 似 的 规 
f& , 80 从 -39.35%o 逐 渐 富 集 到 13.73%0 , BH 从 
-166.8%o 逐 渐 富 集 到 18.43%o ,但 在 整体 上 低 于 剩余 
水 体 中 的 稳定 同位 素 (图 1)。 实 测 大 气 水 汽 中 的 8*0 
在 实验 期 间 的 变化 范围 在 -23.4%o~-13.86%o 之 间 , 平 
均值 为 -18.36%o ; SH 的 变化 范围 在 -154.79%o~ 
-90.11%o 之 间 , 平 均值 为 -122.69%o。 基 于 莱 发 岂 实 
验 模拟 的 大 气 水 汽 中 58"0 的 变化 范围 在 -31.02%o~ 
-7.03‰ 之 间 ,平均 值 为 -18.56% , SH 则 在 
-217.74%o~ -113.22%o 之 间 ,平均 值 为 -153.09%o 。 
从 时 间 序列 来 看 ,模拟 大 气 水 汽 5*0 在 实验 前 期 波 
动 较 小 并 与 直接 测量 得 到 的 大 气 水 汽 3*0 表 现 出 较 
好 的 一 致 性 ,在 实验 末期 模拟 大 气 水 汽 5*0 的 波动 
相对 较 大 。 总 体 来 说 ,在 观测 时 段 内 ,模拟 的 大 气 
水 汽 同位 素 和 实测 的 大 气 水 汽 同位 素 变 化 趋势 基 
本 一 致 。 
2.2 实测 值 与 模拟 值 的 对 比 

根据 蒸发 剩余 水 体 中 的 氧 氧 稳定 同位 素 比 率 ， 
可 以 得 到 8*O FISH 的 线性 关系 , 即 蒸发 线 (MH= 
3.7658"0-28.54, 图 2)。 在 自然 条 件 下 的 蒸发 过 程 
中 ,水体 中 同位 素 的 变化 不 仅 受到 大 气温 度 和 相对 
湿度 的 影响 ,也 受到 大 气 水 汽 中 同位 素 的 影响 ， 
蒸发 线 的 斜率 反映 了 蔡 发 过 程 中 非 平 衡 分 馏 程 度 


的 强 弱 , 斜 率 越 小 表明 非 平衡 分 馏 越 强 。 本 次 实验 
的 蒸发 线 和 斜率 (3.76) 符 合 Cati 提出 的 蒸发 线 斜 率 
通常 在 3.5~6 之 间 的 观点 ,表明 水 体 受 到 了 强烈 的 
非 平衡 分 馏 作用 使 得 蒸发 线 的 斜率 远 低 于 全 球 降 
水 线 的 斜率 。 

相 较 于 实测 的 大 气 水 汽 同位 素 ,模拟 所 得 的 大 
气 水 汽 中 PH 在 整体 上 要 低 于 实际 测量 的 SBH( 图 
Ib) ,因此 模拟 值 分 布 在 实测 值 的 下 方 (图 2)。 整 体 
上 看 ,不 论 实 测 还 是 模拟 的 大 气 水 汽 同 位 素 都 分 布 
在 全 球 大 气 水 线 附近 ,与 Gat 的 观点 基本 一 致 。 
但 是 本 研究 的 实验 地 点 相对 湿度 较 低 (46% ) ,同时 
烈 的 动力 作用 使 得 模拟 和 实测 的 大 气 水 汽 线 的 斜 
率 均 低 于 全 球 大 气 水 线 的 斜率 ,同时 蒸发 亚 水 体 的 
体积 相 较 于 一 般 的 自然 湖泊 等 地 表 水 体 偏 小 ,蒸发 
模拟 的 过 程 更 易 受 到 环境 因子 的 影响 ,使 模拟 得 到 
的 大 气 水 汽 线 的 斜率 更 低 。 

2.3 模型 敏感 性 分 析 

将 相对 湿度 .气温 .蒸发 剩余 比 和 蒸发 剩余 水 
体 同 位 素 组 成 等 参数 进行 敏感 性 分 析 ( 图 3)。 利 用 
改变 参数 后 模拟 得 出 的 大 气 水 汽 同 位 素 与 未 改变 
参数 模拟 得 出 的 大 气 水 汽 同 位 素 的 偏差 来 体现 该 
参数 的 敏感 性 ,偏差 越 大 ,表明 该 参数 对 模型 的 敏 
感性 越 强 。 本 研究 中 ,相对 湿度 的 变化 对 于 模拟 结 
果 的 影响 比较 显著 , 当 相 对 湿度 降低 时 2o(o ,标准 
FE ) 时 ,偏差 为 7.65%o , 而 当 相 对 湿度 增加 2c 时 , 偏 
差 为 3.87%o。 相 对 湿度 降低 产生 的 影响 大 于 相对 湿 
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图 2 2019 年 9 一 11 月 蔡 发 和 由 实验 期 间 模拟 
KAIKA Sae) 实测 大 气 水 汽 (5A) 和 蒸发 
RAK Hs L) HSO MISH 的 关系 
Fig. 2 Correlations between 8"O and 6°H in simulated 


atmospheric water vapor (84.55), measured atmospheric water 
vapor (õa) and remaining pan water (6,) during evaporation 


pan experiment from September to November 2019 


度 增加 时 产生 的 影响 (图 3a) 。 非 平衡 条 件 下 的 动 
力 分 馏 与 相对 湿度 具有 很 高 的 相关 性 ””, 当 相对 
湿度 降低 时 ,动力 作用 加 强 使 得 蒸发 朝 着 更 加 不 平 
衡 的 方向 发 展 ,模拟 的 结果 也 出 现 较 大 的 偏差 。 

在 同等 标准 差 变 化 下 ,温度 对 于 该 模型 的 影响 
相 较 于 相对 湿度 要 小 , 当 温 度 变 化 2 倍 标 准 差 时 , 模 
拟 结果 8*0 偏差 在 3%o 以 内 ,对 于 模拟 准确 性 的 影响 
TEDSEHEE/]NCEEL 3b) ;. YE EL 3c P , RERA EG fI] RTT 
XE CSI , 850 4i 2:5 AY 2E £i BE AI Fi] , EH TP PRIUS AS 
能 超过 1 ,所 以 选择 了 0.3 倍 标准 差 作 为 上 限 ,在 / 增 
加 0.3 倍 标准 差 以 后 ,偏差 增加 的 幅度 也 要 大 于 /7 减 
小 0.3 倍 标准 差 产 生 的 偏差 。 剩 余 水 体 中 同位 素 富 
集 和 贫 化 程度 对 模型 的 影响 相对 比较 对 称 (图 3d)。 
2.4 模拟 可 靠 度 评估 

在 本 次 实验 中 ,对 于 50 而 言 ,RMSE 为 4.5%o， 
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间 模 拟 大 气 水 汽 5"0 在 不 同 相对 湿度 (a) 温度 (b)、 


蒸发 剩余 比 4c) 和 蕉 发 剩余 水 体 j"0(d) 条 件 下 的 敏感 性 分 析 


Fig.3 Sensitivity analysis of simulated 8"O in atmospheric water vapor under various conditions of relative humidity (a), 


temperature (b), remaining ratio (c) and 5"O in remaining pan water (d) during evaporation pan experiment 


from September to November 2019 


MAE 为 3.2%o,MBE X 0.03%0; X} FH 而 言 ,RMSE 为 
38.1%o ,MAE 为 28.3%o ,MBE X 28.3966, Gibson £F” 
在 加 拿 大 北方 的 实验 结果 与 本 实验 相符 ,在 温度 和 
相对 湿度 已 知 的 情况 下 ,利用 蒸发 生 实 验 可 以 得 到 
较为 准确 的 大 气 水 汽 同 位 素 组 成 ,这 在 降水 稀少 的 
地 区 具有 重要 的 价值 。 类 似 的 结论 可 见于 Devi 等” 
和 Vallet-Coulomb 等 的 实验 。 

将 实测 水 汽 8*0 与 模拟 水 汽 8"0 的 差 值 定义 为 
模拟 残 差 ,8380 的 模拟 残 差 值 在 -8%o~16%o 之 间 , òH 
的 模拟 残 差 值 在 0%o~120%o 之 间 。 对 模拟 残 差 与 各 
实验 要 素 进 行 相 关 分 析 发 现 , 温 度 与 38s0 的 模拟 残 
差 之 间 的 相关 性 不 明显 (R=0.08) 且 未 通过 显著 性 
检验 (P>0.05) ,温度 不 是 引起 模拟 残 差 的 主要 原因 
(图 4a)。 而 相对 湿度 与 S8s0 的 模拟 残 差 之 间 的 相关 
性 较 强 (R=0.43) 量 通过 了 显著 性 检验 (P<0.01)， 
两 者 之 间 具 有 负 相 关 性 (图 4b)。 相 对 湿度 在 40% 
以 上 时 ,580 的 模拟 残 差 数据 点 分 布 较为 集中 ,而 在 
相对 湿度 小 于 40% 时 模拟 残 差 的 数据 点 分 布 则 更 
加 离散 ,880 模拟 残 差 的 极 大 值 也 出 现在 相对 湿度 


小 于 40% 的 这 一 区 间 内 , 且 对 应 了 相对 湿度 的 最 小 
值 。 在 低 相 对 湿度 的 情况 下 ,水 体 在 蒸发 过 程 中 会 
经 历 更 强烈 的 非 平 衡 分 馏 过 程 ,使 得 模拟 结果 出 现 
较 大 的 偏差 。 水 汽 压 与 5*0 模 拟 残 差 呈 负 相关 , 相 
关 性 低 于 相对 湿度 和 8*0 模 拟 残 差 的 相关 性 ,数据 
点 的 分 布 也 较为 分 散 (图 4e)。f 与 5*0 模 拟 残 差 呈 
HAASE (020.39) ,在 f<0.3 以 后 ,5"*0 的 模拟 残 差 值 
分 布 较为 离散 且 出 现 了 极 大 值 和 极 小 值 ,而 在 0.3 < 
了 <0.8 的 部 分 ,8*0 模 拟 残 差 的 分 布 比较 集中 ,多 在 
0%o 附 近 ( 图 4c)。78 较 小 时 表明 蒸发 剩余 水 体 较 少 ， 
易 受 到 水 体 和 大 气 之 间 同 位 素 交换 作用 的 影响 , 进 

影响 蒸发 严实 验 的 结果 。 蒸 发 剩余 水 体 880 与 
8^0 模拟 残 差 具有 一 定 的 相关 性 ( 尼 =0.37) ,在 实验 
前 期 ,蒸发 剩余 水 体 中 同位 素 的 升 高 ,模拟 残 差 从 
负 值 逐渐 过 渡 到 正 值 , 但 总 体 仍 比较 接近 0。 随 着 
蒸发 剩余 水 体 880 大 于 8%o 以 后 ,模拟 残 差 的 分 布 变 
得 更 加 离散 (图 4d)。 在 莱 发 实验 的 后 期 ,蒸发 了 中 
水 量 减少 后 更 易 受 到 周 于 环境 变化 的 影响 ,使 得 薰 
发 过 程 变 得 不 稳定 。 同 时 ,在 蒸发 后 期 ,蒸发 由 中 的 
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Fig. 4 Correlations between residual of simulated õ'*O in atmospheric water vapor and temperature (a), relative humidity (b), 


水 体 中 重 同 位 素 已 经 变 得 相对 富 集 ,实验 后 期 8s0 
的 模拟 残 差 较 大 且 较 为 分 散 。 和 气压 与 Ss0 的 模拟 
残 差 之 间 则 没有 表现 出 明显 的 相关 性 (R=0.001) 
(图 4f)。 
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remaining ratio (c), ^O in remaining pan water (d), water vapor pressure (e) and atmospheric pressure (f) 


during evaporation pan experiment from September to November 2019 
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8 于 的 模拟 残 差 与 各 要 素 间 的 相关 关系 如 图 5 
Bron ,整体 上 与 图 4 中 88s0 的 结果 类 似 , 但 8H 的 模 
拟 残 差 在 数值 上 要 大 于 380 的 模拟 残 差 。 此 外 ,5°H 


的 模拟 残 差 与 湿度 (R=0.51) ,水 汽 压 (R=0.53) 28 
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图 5 2019 年 9 一 11 H 2& A MSC ds T 
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Fig. 5 Correlations between redidual of simulated ôO in atmospheric water vapor and temperature (a), relative humidity (b), 


evaporation pan experiment from September to November 2019 


remaining ratio (c), &'*O in remaining pan water (d), water vapor pressure (e) and atmospheric pressure (f) during 
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相关 性 要 强 于 8*0 的 模拟 残 差 与 上 述 要 素 间 的 相关 
性 。 在 湿度 低 于 40% 和 < 0.3 以 及 莹 发 剩余 水 体 中 
同位 素 富 集 以 后 ,SH 的 模拟 残 差 均 出 现 了 较 大 的 
离散 ,但 在 相对 湿度 高 于 40% 和 0.3<f<0.8 以 及 燕 
发 实验 前 期 ,YH 的 模拟 残 差 与 上 述 三 要 素 则 体现 
出 较 好 的 相关 性 。 总 体 来 说 ,相对 湿度 KERR 
比 和 蒜 发 剩余 水 体 中 的 氧 氧 同位 素 三 者 都 与 3*0 和 
8 五 的 模拟 残 差 存在 明显 的 联系 。 


3 结论 


本 研究 基于 同位 素质 量 平衡 和 Craig-Gordon 线 
性 阻力 模型 设计 了 莹 发 亚 蒸 发 实验 ,模拟 了 与 茸 发 
过 程 同期 的 大 气 水 汽 同位 素 组 成 ,结合 实际 监测 的 
大 气 水 汽 同 位 素数 据 , 对 该 模型 模拟 的 可 靠 性 做 了 
评估 ,得 出 了 以 下 结论 : 

(1) 在 实验 期 间 内 ,基于 蒸发 刀 实 验 模拟 的 大 
气 水 汽 中 8*0 与 实际 测量 的 大 气 水 汽 中 5*0 相 比 ， 
RMSE 为 4.49%o ,MAE 为 3.23%o, MBE 为 0.03%o , 模 
拟 得 到 的 大 气 水 汽 同位 素 变化 可 以 反映 出 大 气 水 
汽 同 位 素 的 变化 。 

(2) 通 过 对 模型 的 敏感 性 分 析 发 现 , 在 同等 标 
准 差 变化 下 ,温度 对 于 该 模型 的 影响 相 较 于 相对 湿 
度 要 小 。 蒸 发 剩余 比 减 小 和 增 大 时 ,3580 偏差 的 变 
化 幅度 不 同 ,蒸发 剩余 比 增 大 时 结果 更 加 敏感 。 

(3) 相 对 湿度 、 蒸 发 镜 余 比 和 剩余 水 体 的 同位 
素 值 都 与 水 汽 5s0 ASH 的 模拟 残 差 存在 联系 。 当 
相对 湿度 低 于 40% 时 ,模拟 残 差 的 分 布 更 加 离散 。 
在 实验 末期 蒸发 剩余 比 减 小 以 后 ,模拟 残 差 相 较 于 
实验 前 期 而 言 更 大 , 较 少 的 水 量 和 蒸发 剩余 水 体 中 
相对 富 集 的 氨 氧 稳定 同位 素 对 应 着 较 大 的 水 汽 同 
位 素 模拟 残 差 。 
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Simulation of stable hydrogen and oxygen isotopes in atmospheric water vapor 


based on an evaporation pan experiment 


LEI Shijun’, WANG Shengjie", ZHU Xiaofan, ZHANG Mingjun" 
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Abstract: The stable hydrogen and oxygen isotopes in atmospheric water vapor reflect the key processes of water 
transport, mixing, and phase change in the atmosphere. The stable isotopic compositions of hydrogen and oxygen 
in atmospheric water vapor can be simulated using an evaporation pan experiment, but a measurement-based assessment 
of model reliability is necessary. In this study, an evaporation pan experiment was conducted in Lanzhou, Gansu, 
from September to November 2019. The daily atmospheric water vapor isotopes were simulated using the Craig- 
Gordon linear resistance model and the water isotopic mass balance method in a natural evaporative state. Online 
observations of atmospheric water vapor isotopes were performed using a water vapor isotope analyzer and then 
the pan-based simulations were verified. The simulated O in atmospheric water vapor and the directly measured 6 
8O were in good agreement with a root mean square error of 4.5%, a mean absolute error of 3.2%o, and a mean 
bias error of 0.03%o .The residual of simulated O values were correlated with relative humidity (R’=0.43), 
remaining ratio (R^—0.39), and the isotopic values of the remaining water bodies (R =0.39). The model was highly 
sensitive to relative humidity: lower relative humidity conditionstypically corresponded to a higher error. Under 
the conditions of low remaining ratio and isotopically enriched remaining water at the end of the evaporation 
experiment, the simulation residual was much larger than that at the beginning of the experiment. 


Keywords: stable hydrogen and oxygen isotopes; evaporation pan; simulation; atmospheric water vapor 


